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The Sommerfeld and Brillouin forerunners generated in a single-resonance absorbing medium 
by an incident step-modulated pulse are theoretically considered in the double limit where the 
susceptibility of the medium is weak and the resonance is narrow. Combining direct Laplace- 
Fourier integration and calculations by the saddle-point method, we establish an explicit analytical 
expression of the transmitted field valid at any time, even when the two forerunners significantly 
overlap. We examine how their complete overlapping, occurring for shorter propagation distances, 
originates the formation of the unique transient currently named resonant precursor or dynamical 
beat. We obtain an expression of this transient identical to that usually derived within the slowly 
varying envelope in spite of the initial discontinuity of the incident field envelope. The dynamical 
beats and Ott pulses generated by ultra-short incident pulses are also briefly examined. 
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I. INTRODUCTION 

About one century ago, in order to remove an appar- 
ent inconsistency between classical theory of waves and 
special relativity, Sommerfeld and Brillouin studied in de- 
tail the linear propagation of step-modulated light pulses 
in an absorbing medium with a single absorption line 
P, i|. They found that, at propagation distances such 
that the medium is opaque in a broad spectral range, 
the transmitted field consists of two successive transients 
preceding the establishment of the steady-state field at 
the carrier frequency ujc of the pulse (the "main field"). 
They naturally named the transients "forerunners". To 
be definite, we will conserve this name in the following, 
reserving the name of precursor to the unique transient 
occurring when the two forerunners completely overlap. 
Sommerfeld and Brillouin showed that the first forerun- 
ner (the "Sommerfeld forerunner") and the second one 
(the "Brillouin forerunner") respectively involve frequen- 
cies large and small compared to the frequency ujq of 
the absorption line (that is in the spectral regions where 
the medium has some transparency). They proved that 
the front of the first one propagates at the velocity c of 
light in vacuum, in agreement with special relativity, and 
pointed out that the second one approximately moves at 
the group velocity at zero frequency. The shape of the 
forerunners was derived by means of classical complex 
analysis fl'| and the newly developed saddle point method 
[2]. Later Brillouin used the stationary phase method to 
study his eponym forerunner when its formation is domi- 
nated by dispersion effects Q . Following these pioneering 
works [J| , the forerunners became a canonical problem in 
physics and entered reference textbooks in electromag- 
netism @|. More rigorous solutions, correcting the re- 
sults obtained by Sommerfeld and Brillouin, were derived 



'Electronic address: |bernard . segard@univ-lille- 1 . fr | 



by means of uniform asymptotic methods. See, e.g., 0- 
Il4j . The problem was also studied by a purely t emp oral 
approach instead of the usual spectral approach [I3 ■ An 
abundant bibliography on forerunners can be found in 
[T^ . For most recent studies related to the forerunners 
in the sense of Sommerfeld and Brillouin, see 

The expressions of the forerunners obtained in the 
most general case [lo| are tremendously complicated 
and not explicit. A good insight on the physics of fore- 
runners is fortunately obtained by examining particular 
cases. We recently established simple analytical forms 
for the Sommerfeld and Brillouin forerunners when the 
propagation distance is such that the two forerunners are 
far away from each other [20]. In the present paper we 
study the problem when the electric susceptibility of the 
medium is weak and the resonance is narrow. This dou- 
ble condition is generally satisfied for dilute media in the 
entire optical domain but also for dense media in the X 
and 7 spectral regions. We show that it is then possible to 
obtain explicit analytical expressions for the forerunners, 
valid even when the two forerunners significantly overlap, 
and to study how their complete overlapping originates 
the precursors actually observed in optics |2l| - |27| or the 
dynamical beats evidenced in the experiments of nuclear 
coherent forward scattering [28, 29J. The arrangement of 
our paper is as follows. In Sec. [Ill we give the transfer 
function of the medium in the considered limit of weak 
susceptibility and narrow resonance. We establish in Sec. 
mil the expressions of the Sommerfeld and Brillouin fore- 
runners. We study in Sec. IIVI the evolution of the fore- 
runners towards a unique transient (optical precursor or 
dynamical beat) and show that the expression of the lat- 
ter is identical to that obtained within the slowly vary- 
ing envelope approximation even when the envelope of 
the incident field is initially discontinuous. Finally the 
dynamical beats and Ott pulses generated by ultra-short 
incident pulses are briefly revisited. We conclude in Sec. 
|V]by summarizing our main results. 
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II. TRANSFER FUNCTION OF THE MEDIUM 

We consider a one-dimensional electromagnetic wave 
propagating in the ^-direction through an isotropic and 
homogeneous medium. Its electric field is assumed to be 
polarized in the x-direction {x,y,z : Cartesian coordi- 
nates). We denote e{0,t) the field at time t for z = 
(inside the medium) and e{z,t) its value after a propa- 
gation distance z through the medium. In a spectral ap- 
proach the medium is characterized by its transfer func- 
tion H(z,uj) relating the Fourier transform E{z,u;) of 
e(z, t) to that £^(0, uj) of e(0, t) 



E{z,uj) = H{z,uj)E{0,uj). 



(1) 



The transmitted field, inverse Fourier transform of 
E(z, w), reads as: 



e{z,t)= I H{z,u)EiO,u)e"^'^. (2) 



In all the following, we take for t a retarded time equal 
to the real time minus the luminal propagation time z/c 
(retarded-time picture). We have then 

i7(z,cj) ==exp|-2^:!^ [v/TTxH-i]}- (3) 

Here x(a;) is the complex electric susceptibility of the 
medium at the frequency uj and -^/l -I- xi^) is its complex 
refractive index. x('^) being given, Eq.Q can be numer- 
ically solved by fast Fourier transform (FFT) but has no 
exact analytical expression, even for the simplest forms 
of x{^)- Following Sommerfeld and Brillouin, most au- 
thors have considered a Lorentz medium consisting of an 
ensemble of damped harmonic oscillators with the same 
resonance frequency uq and the same damping rate 7. 
Its susceptibility reads as: 



x(w) 



2iju} 



(4) 



where ujp is the so-called plasma frequency whose square 
is proportional to the number density of absorbers. A 
very similar expression of the susceptibility is obtained 
for the two-level medium usually considered in quantum 
optics Without doing the usual rotating wave and 
slowly varying envelope approximations, we get: 



■ uiq — 2i^uj — 7^ 



(5) 



shown in particular that the very first beginning of the 
transmitted field always propagates without distortion at 
the velocity c and that its total area (to distinguish from 
that of its envelope) is conserved during the propagation. 

The problem of the forerunners is greatly simplified in 
the double limit considered in the present paper where 
7 ^ (narrow resonance) and |x(w)| ^ 1 at every fre- 
quency (weak susceptibility). Due to the first condition 
the susceptibility of the two-level medium equals that of 
the Lorentz medium (7^ negligible with regard to ujq). 
The second condition is fulfilled when ojp <C 27^0, that 
is, owing to the first condition, when u)p <^ ljq. The 
absorption coefficient for the amplitude aicj) then takes 
the simple form a{uj) ~ — ^Im[x(a;)]. It is everywhere 
small compared to the wavenumber fc(w) — uj/c and is 
maximum for iv « iojQ with a(±aJo) = ao ~ ujp/{A'yc). 
It is convenient to characterize the propagation distance 
by the corresponding optical thickness on resonance aoz. 
The transfer function then reads as: 



iJ(z, uj) w exp 



iujzx{(^) 
2^ 



exp 



2iaozjui 



2i7w 



(6) 

Eq.® is obtained by expanding a/1 -|- xi'^) the first 
order in x(w). Its upper limit of validity is given by the 
condition : 



LJZ 
C 



< 1 



(7) 



This condition has to be fulfilled at the frequencies u for 
which the medium has some transparency, say for which 
|i?(cj)| > 10^'*. For the parameters considered hereafter, 
numerical simulations show that Eq.Q is over-satisfied 
for resonance optical thickness uqz up to 10^. 

An alternative form of H{z,uj) can be obtained by ex- 
panding the exponent in Eq.([n]) in partial fractions. It 
reads as H{z,uj) ~ H^{z,uj) ■ H^{z,uj) where 



H±{z, uj) = exp 



laozj 



1 ± ij/ujQ 



(8) 



Under this form, the transfer function is very similar to 
that encountered in the study of optical precursors in a 
medium with a transparency window between two ab- 
sorption lines [3^. There is obviously some analogy be- 
tween the two problems. It should be noticed, however, 
that the transfer function considered in [s^l was associ- 
ated with the envelope of the electric field whereas that 
given by Eq.® and Eq.® is associated with the field 
itself. 



with ujp — 2ujop /hsQ. In this expression p is the dif- 
ference of population per volume unit between the two 
levels at thermal equilibrium, /i is the dipole moment ma- 
trix element of the transition and 7 is the relaxation rate 
for the coherence. General properties of the transmitted 
field resulting from Eqs.(I3][S]) are discussed in [20J. It is 



III. SOMMERFELD AND BRILLOUIN 
FORERUNNERS 

As Sommerfeld and Brillouin and most authors, we 
consider in this section a step-modulated incident field 
of the form e{0,t) = sin (wct) • unit) where unit) is the 
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Heaviside unit step function and Wc is the carrier fre- 
quency. Eq.® can then be reduced to 



e(z, t) — Im 



(9) 



where F is a straight hne parallel to the real axis passing 
under the pole at w = ojc- For the large propagation 
distances at which the forerunners are discernible, the 
medium is opaque in a broad spectral region and the 
ranges of action of i?+ and H- overlap. We assume here 
that Wc lies in the opacity region. The transmitted field 
e(z, t) then only contains high (w > uj^ and low (w < Wc) 
frequencies, respectively associated with the Sommerfeld 
and Brillouin forerunners. We write it : 



where J„(s) designates the Bessel function of the first 
kind of index n. For 2ujQy/tBt ^ 1 the Bessel function 
may be replaced by its asymptotic form [s^l to yield 



ei{z,t) « — 

UJQ 



2t cos(2a;o-y 



(15) 

Ji(s) equaling its asymptotic form for s w 2.47, the fields 
given by Eg. p^ and Ea. p4| are equal at the time t — 
h ~ 1.52/ {uj^tB). 

Similarly we get for the beginning of the Brillouin fore- 
runner [35| 



e2iz,t) « AAi(-60e-2Tt73. 



(16) 



e{z,t) = ei{z,t) + e2{z,t) 



(10) 



where ei{z,t) and e2{z,t) respectively stand for the first 
(Sommerfeld) forerunner and the second (Brillouin) fore- 
runner. The forerunners are determined both by the 
frequency-dependence of the absorption of the medium 
and by its dispersion. The latter may be characterized 
by the group delay Tg{z,uj) = —d^/duj where $(z,a;) 
is the argument of H{z,lo). For the high and low fre- 
quencies associated with the Sommerfeld and Brillouin 
forerunners, respectively, we get the asymptotic forms 



Tg (z, |a;| > Wo) 



2ao^7 



(11) 



and 



Tg (z, |w| < Wo) ~ is + 



tB 1 



where 



tB = 



2aQZ"f 



= Tg{z,0)-Tg{z,(X)) 



(13) 



tg is obviously indicative of the time-delay of the Bril- 
louin forerunner with respect to the Sommerfeld forerun- 
ner and provides a good time scale for the study of both 
forerunners. 

Eq. lfTTj) and Eq.(fT2|l show that the Sommerfeld fore- 
runner will start at the retarded time t = with a in- 
finitely large instantaneous frequency whereas the latter 
is vanishing for the Brillouin forerunner at t « i^. As 
already numerically evidenced (see Fig. 9 in [13), the be- 
ginning of the forerunners will thus be well reproduced 
by using asymptotic forms for H{z,uj) and E(0,uj). The 
transmitted field can then be calculated by direct inte- 
gration of Eq.([9]) by means of standard Laplace- Fourier 
procedures [20]. 

For the beginning of the Sommerfeld forerunner, we 
get [H: 



ei(z,t) 



•e-^^'uH{t) (14) 



Here b = [uj^/ (Gao^T)] = KV (S^s)] t' = t - tg 
and Ai(s) designates the Airy function. This solution 
is physically acceptable if and only if e2(z,t) ~ for 
t < 0. This is achieved when Ai{btB) ~ 0, that is when 
aQZ ^ W0/7. On the other hand, for bt' ^ 1, we get 
from the asymptotic form of Ai(— s) (s^ 



,1/3 



62(2, t) 



,1/2 



sm 



2wo^'3/2 

3(3ts)'/' 



(17) 

Ai(— s) equaling its asymptotic form for s w 1.42, the 
fields given by Ea. (|17p and Eq. p^ are equal at the re- 
tarded time i = ^2 ~ is + 1.42/5 w is 2.05 (tjs/wg)^^^. 

The previous expressions of the forerunners are only 
valid up to a finite time, the longer the larger the prop- 
agation distance [i^l- We find that the corrections to 
Eq. (fT4|) [Eq.(|16p] due to the next terms in the asymp- 
totic expansion of ln[i7(z,w)] and of E{0,u>) in powers 
of 1/cj [w] are negligible up to i = ti \t = t2] when the 
condition a^z ^ wo/7 is satisfied. This condition thus 
suffices for the vahdity of Eq.dlll) [Eq.([Tl])] for < t < ii 
[0 < t < t2]. 

For t > ti {t > t2) , the Sommerfeld (Brillouin) fore- 
runner is calculated by means of the basic saddle point 
method as used by Brillouin [2, 4]. Introducing the phase 
function 5'(z,w) = iutt + In [H{z,uj)], Eq.® is rewritten 
as 



e(z, t) — Im 



exp [4' (z, uj)] duj 



2iTT 



(18) 



The integral is calculated by deforming the straight line 
r in a contour travelling along lines of steepest descent 
of the function 5'(z,a;) from the saddle points where 
d'^ /duj = [13 . The contribution of a nondegenerate 
saddle point at 0;^ to the integral reads as 



a{ujs) 



exp (ws) + it 



i {ujs ~ Wc) v/27r*" (w,) 



(19) 



where is the angle of the direction of steepest descent 
with the real axis and ^E*" (ws) is a shortcut for jduP' 
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at = a;^. We may disregard the presence of the pole 
at ujc- The corresponding residue is indeed neghgible 
when, as assumed here, lOc hes in the opacity region of 
the medium. The determination of the saddle points is 
very simple in the double limit considered in the present 
paper. Due to the weak susceptibility hypothesis, the 
equation giving the complex frequencies of the saddle 
points is only of fourth degree (instead of eighth degree 
in the general case) and the narrow-resonance condition 
allows us to solve this equation at the lowest order in 7. 
Regrouping the four solutions in two pairs, we easily get 



where 



= ±a;„ + i7„ with n = 1,2 and 



: UJ0\ 1 + 



tB_ 

2t 



(-1)" 



8t 



1 + f 

IB 



(20) 



(21) 



(22) 



The pair associated with wi (high frequency) obviously 
originates the Sommerfeld forerunner and yields 



ei(z,t) = Im [fli (tj+) + fli (u^ )] 



(23) 



ei{z,t) 
I 

X 



[ujit + wiw^is/ {ujI - ^2) - 37r/4] 



\-3 



where 



-71* 



(24) 



(25) 



For t tB, Eq.lgl]) is reduced to Eq.lUn]), which itself 
fits Eq. (|14p at the time t ^ ti. This means that the 
fields given by Eq. ([24)) and Eq. ([T4|) will fit together in ti 
if ti ^ ts. It is easily shown that this is achieved when 
a^z wo/7- The combination of Ea. ([H)) for t < ti and 
Ea. (|24p for t > ti then provides an analytical expression 
of the Sommerfeld forerunner valid at every time. 

A similar calculation for the pair of saddle points asso- 
ciated with UJ2 (low frequency) yields the Brillouin fore- 
runner 



62(2:, t) 
/ 

X 



UJoUJc 



TT \0J, 



sin [iJ2t - ijJ2i^ltBl (wq - ^^2) + '''/4] 



-tB 



(^2+1:^^0) +{uJ2-(^oy 



-72*' 



(26) 



1 



^-3 

tB 



(27) 



For t' ^ t - tB <^ tB, Eq.(l26l) is reduced to Eq.(fT7|. 
which itself fits Eq. p^ at the time t2. lit2 — tB^tB (a 
condition also satisfied when a^z ^ ujq/j), the Brillouin 
forerunner will thus be well reproduced by Eq. (|16p for 
t < t2 and by Eq.dH]) for t > t2. 

To summarize, when the propagation distance is such 
that a{ujc)z ^ 1 (opacity condition for the main field) 
and that a^z 3> wo/7, the transmitted field is simply 
the sum of the Sommerfeld and Brillouin forerunners for 
which we have obtained piecewise analytical expressions 
valid in the entirety of the time domain where they have 
significant amplitude. 
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Figure 1: (Color online) Transmitted field e{z,t) for an in- 
cident field e{Q,t) = sm{uict) ■ Uu{t)- Parameters 7 = ojp = 
ajo/200, LOc ~ 100 and aoz — 10^, for which aoZ"f/ujo — 500, 
iootB = 1000 and 'yts ~ 5 . The full line (dashed line) is the 
exact numerical solution (our approximate analytical solu- 
tion). For the considered parameters, the first (Sommerfeld) 
and second (Brillouin) forerunners are well separated. Up- 
per (lower) inset: enlargement of the beginning (tail) of the 
Sommerfeld forerunner. 

To illustrate our results, we consider first the case 
where the two forerunners are well separated. In order to 
meet the weak-susceptibility and narrow-resonance con- 
ditions, we take j = ujp — wo/200. Its damping time 
being of the order of I/7, the Sommerfeld forerunner will 
not overlap the Brillouin forerunner if ^ I/7, that is 
when aoz ^ cjq/ (27^). Figure [1] shows the result ob- 
tained when "ftB = 5, a value attained at a propagation 
distance z such that aoz = 10^ for which wgis = 1000, 
ujoti fa 1.52 X 10^"^ and ujq {t2 ~ tB) ~ 20.5. As expected, 
our piecewise analytical solution perfectly fits at every 
time the exact numerical solution obtained by using the 
transfer function given by Eq.([3]) and Eq.Q without any 
approximation. It is worth noticing that the maximum 
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of the Brillouin forerunner occurs at a time shorter than 
^2 and that 7 {t2 — is) ^ 1. Eq. (|16P then shows that the 
corresponding amphtude is proportional to & cx z"^/'^ 
and inversely proportional to iOc, no matter the location 
of iOc in the opacity region, inside or outside the anoma- 
lous dispersion region. We also remark that, for the large 
optical thickness considered in this example, the asymp- 
totic form of the Sommcrfeld forerunner given by Eq. (|14p 
holds much beyond ti. This equation even provides a 
good estimate of the maximum amplitude of the forerun- 
ner that occurs at t « 1/ (87) and is proportional to Wc- 
We have checked all these points by comparing the ex- 
act numerical results obtained for lJc equal to wq (Fig[T]), 
LU0V2 and ujo/^/2. 
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Retarded Time (units of I/cOq) 



Figure 2: (Color online) Same as Fig. [T] for aoz = 1000 
(aoz"//ijJo — 5 , uJotB = 10 , and jtB = 1/20). The two 
forerunners significantly overlap but remain discernible. The 
upper inset clearly shows the oscillations of the Sommerfeld 
forerunner superimposed over the slow rise of the Brillouin 
forerunner. In the far wing (lower inset) appear clean beats 
that anticipate what will become the optical precursor. 

In the previous example, the condition of validity of 
our analytical results {agz ^ ojq/"/ ) was over-satisfied. 
Figure [5] shows, other things being equal, the result 
obtained for a propagation distance 100 times shorter, 
for which aozj/uo — 5, uJots — 10, ujoti ~ 0.15 and 
<^o(i2 — is) ~ 4.4. Though the condition of validity 
of our approximations is then marginally satisfied, our 
piecewise analytical solution continues to fit very well the 
exact numerical solution. Since is is now much shorter 
than 1/7 {jtB = 1/20), the Sommerfeld forerunner sig- 
nificantly overlaps that of Brillouin. It remains quite vis- 
ible for t < ts where it is superimposed to the slow rise 
of the latter (see upper inset of Figl5]). For t > tB , beat- 
ings between the two forerunners are observed [Ol • They 
are well developed when the instantaneous frequencies 
uji and UJ2 of the forerunners are close and their ampli- 
tudes are comparable. This occurs when y/t/ts ^ 1 [see 
Eqs.il^ni HH HH)]. In the present case, the corresponding 



times are long compared to the damping time I/7 and 
the amplitude of the beats is weak (see lower inset of 
FigEl). 
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Figure 3: (Color online) Same as Fig. [T] for aoz — 200 
(0:027/^0 = 1 and LdotB = 2). A sole oscillation of the Som- 
merfeld forerunner is visible and clean beats appear sooner 
with significant amplitude. The analytical solution (dashed 
line) is only given in the time domain t > t2 where it is re- 
duced to the saddle point solution (see text). 

For smaller propagation distances (that is for shorter 
times ts), clean beatings are expected to appear sooner 
with larger amplitudes. This is illustrated Fig[2] obtained 
for aQZ"f/ujQ = 1, that is for aoz — 200 with our pa- 
rameters. We have then looIb = 2, woti ~ 0.76 and 
ujQt2 ~ 4.58. For t < t2, a single oscillation of the Som- 
merfeld forerunner is still visible. Since the condition 
00^7/1^0 3> 1 is not satisfied, it should be noticed that 
our piecewise solutions for ei(z, t) and 62(2, t) fail in this 
time domain. On the other hand, the latter continue to 
perfectly fit the exact numerical solutions for t > t2, & 
domain where they are reduced to the solutions given by 
the saddle point method [Eqs. d^ B^ ]. As we shall see 
later, these solutions remain a good approximation for 
t > i2 as long as the medium is opaque at the carrier 
frequency. 

Finally when the propagation distance still decreases 
to become such that aoz^/ujQ <C 1 and thus wq^b ^ 
1, the two forerunners completely overlap to originate 
the unique transient currently called resonant precursor 
[itI [37|. Sommerfeld-Brillouin precursor [l^ or dynami- 
cal beat f29'|. It is then possible to obtain an analytical 
expression of the transmitted field valid at every time, no 
matter the optical thickness a^z. 



IV. OPTICAL PRECURSORS OR DYNAMICAL 
BEATS 

This study is facilitated by writing the field under the 
form e(z, t) = Re [eo(2:, t) e"^"*] where eo(z, t) is the com- 
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plex envelope of the field in a frame rotating at the angu- 
lar velocity -l-Wo- We consider an incident field e{0,t) = 
Re [f{t) e^'^"*] that generalizes that considered in Sec lIIII 
Its complex envelope reads as cq (0,t) = f{t) e*^* where 
A = ujc — ujQ is assumed to be small compared to luq. 
Eq.© then leads to 



eo{z,t) 



(28) 



where F{Ld) is the Fourier transform of f{t). When 
00^7/1^0 <C 1 as considered in all this section, the partial 
transfer functions H+{z,uj) and iJ_(z,w) defined Eq.® 
significantly differ from 1 only in narrow domains around 
+iOo and — wq, respectively. It is then justified to make 
the so-called rotating wave approximation [3l!| and to ap- 
proximate H{z,Lu) in Ea. (P5| by H+{z,uj). Translating 
the frequencies by —loq in the integral, we get 



eoiz,t) = / Ho{z,Lu)F{Lu-A)e'' 



> duj 
2^ 



(29) 



where 



Hoiz,uj) = exp 



-aozj 



1 



7-1-10; / 
^ exp 



"f + iuj 



(30) 



Here Hq{z,uj) is nothing else that the transfer function 
for the field envelope and characterizes the medium in- 
dependently of the incident field. The corresponding im- 
pulse response is easily obtained by inverse Laplace trans- 
form and reads as 



ha{z,t)^5{t) 



t 



Ji{2y/^^;i^e-'^'uH{t) (31) 



where 6{t) is the Dirac delta function. On the other hand 
F(aj— A) is the Laplace-Fourier transform of /(t)e*'^* and 
we get from Eq. (|29|) : 



eo{z,t)^ho{z,t)(E) [f{t) 



JAtl 



(32) 



Combined with the relation e{z, t) — Re [eo{z, t) e''^°*] , 
Eq. (|32p enable us to determine the transmitted field 
for arbitrary modulation of the envelope of the incident 
field. When the latter is step modulated, that is when 
f{t) cx Unit), the convolution of Eq. (|32p is easily calcu- 
lated. Replacing this result in the expression of e{z,t), 
we finally obtain the transmitted fields for the incident 
fields sin (ujct) Unit) and cos (ujct) unit). They respec- 
tively read as Im [e (z, t)] and Re [e{z, t)] with 



In agreement with the Feynmann analysis of the ab- 
sorption and dispersion phenomena in linear media (38| . 
Eq. ((33|) makes explicit in both cases that the transmit- 
ted wave is the sum of the incident wave as it would 
propagate in vacuum and of the secondary wave radi- 
ated by the polarization induced in the medium, initially 
of zero amplitude 0, [so']. Results equivalent or ana- 
log to those given En s. ([311(55)) were established in the 
past [m, Hll, fsol - lill lisHi^ However they were generally 
obtained in the frame of the slowly varying envelope ap- 
proximation (SVEA) [3l|. As soundly remarked in [iq . 
its use is quite disputable when the envelope is initially 
discontinuous. SVEA is not made in our calculations. 
In order to check the validity of the latter, we compare 
FiglHthe fields derived from Eq. p3|) to the exact numeri- 
cal solution for aoz — 20, with 7 and uip as in Sec lIIII We 
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Figure 4: (Color online) Transmitted field (optical precur- 
sor) obtained with incident fields (a) cos(ajct) • unit) and (b) 
sm{ujct) ■ UH{t) for a^z = 20 (0027/0^0 = 1/10 ). In both 
cases, the analytical solution obtained within the rotating 
wave approximation (dashed line) very satisfactorily fits the 
exact numerical solution (full line). 

have then a{)Z^ /loq = 1/10 . The agreement between the 
two solutions is very satisfactory, even when the incident 
field itself is discontinuous (curve a). The discrepancy is 
everywhere smaller than 00^7/^0, which is the order of 
magnitude of the deviation of H±{z, T^o) from unity. 

Although precursors are not generated in this case, we 
incidentally mention that Eq. p3p admits an explicit solu- 
tion when the incident field is significantly detuned from 
resonance (A 3> 7,7ao2;). The exponential e"^* is then 
rapidly variable compared to the rest of the integrand 
and a simple integration per part yields 



e{z,t) 





laaz'j ( 







X e'"=*Ui/(i) (33) 



e{z,t) ? 
002:7 



sin (uJct) 
cos {uJct) ' 



Ji (2^0:027^) 
V<^ozjt 



-It 



cos (woi) 



(34) 
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where the multiphcation by unit) has been omitted for 
simphcity. The first term in Ea. (|34p is the incident field 
which is transmitted with neghgible attenuation when 
A ^ jaQz) whereas the second term evidences a beat 
between the incident field and the field reemitted by 
the medium at its eigenfrequency. Since Ji (s) —> s/2 
for s — > 0, the initial amplitude of the beat is zero, as 
expected. A convincing experimental demonstration of 
such beats can be found in [49|. When aoz <C f , the beat 
is reduced to [cos (uJct) — e~'^* cos (ujot)]. This con- 
dition is approximately met in the experiments reported 
in dOl. 

Precursors are obtained in the opposite case where 
the incident field is resonant or quasi resonant (A ^ 
7,7aoz). We restrict the analysis to this case in the fol- 
lowing. The envelope eo{z,t) of the transmitted field is 
then real and can be written as 



ft 



eo{z,t) 



1 — agz 



Ji {2y/^^) 



aozb 



e-^de unit). 



(35) 

It takes a simplified form when the Bessel function in 
the integral evolves rapidly with respect to the expo- 
nential, that is when a^z ^ I, A/7. Again by an 
integration per parts, Eq. (|35p then yields eQ{z,t) ~ 
Jo {2^aoz^t) e~^*UH{t) and the transmitted field for 
e(0, t) — sin (wct) unit) reads as: 



e(z, t) Jo (2-v/ao^7^) sin (wot) e '^*UH{t)- 



(36) 



It consists of successive lobes of decreasing amplitude 
and increasing duration, separated by zeroes of amplitude 
occurring at times t = Jq^/ {Aaozj), where j„p is the 
zero of order p of J„ (s). Fig. [5] (curve a) shows that 
this approximate analytical solution satisfactorily fits the 
exact numerical solution. As announced in Sec lIIII and 
shown FiglS] (curve b), this is also true for t > t2 for the 
saddle point solution. This simply results from the fact 
that ^2 3> with ujots 1 and that the two solutions 
have then the same asymptotic form for t ^ ts- 

When aoz decreases, our saddle point approximation 
becomes worst and worst owing to the coalescence of wi , 
0J2 and ujQ in a time domain where exp (— 7^) = 0(1) 
whereas the rotating wave approximation becomes better 
and better. Correlatively the number of lobes in e{z,t) 
and of zeroes for the amplitude decreases. Of special in- 
terest is the case where there is a single zero of amplitude. 
The integral in Eg. ([55]) having its first maximum (abso- 
lute maximum) for 2^/aozjt = jn, this will obviously 
occur when this maximum is equal to 1, that is when 



/ 



Ji (2Vd 



exp 



aoz 



dO = 1. 



(37) 



This equation in a^z is easily solved by numerical proce- 
dures to yield a^z w 2.80 , the zero of amplitude being 
attained at the time t — tc such that ^tc = jii/ {Aa.i^z). 
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Figure 5: (Color online) Extended view of the precursor ob- 
tained in the conditions of Fig[4]with e(0, t) — sin(cjct) ■ UH{t) ■ 
The exact numerical solution (full line) is compared (a) to the 
simple solution given by Eq. (|36p and (b) to the saddle point 
solution for t > t2 (dotted line). 



This solution is general and does not depend on a particu- 
lar choice of parameters. Figure [5] shows the transmitted 




12 3 

Retarded Time (units of l/y) 

Figure 6: Transmitted field as a function of 74 for aoz — 2.80. 
The precursor has a sole lobe and its amplitude falls to 
for = ^tc ~ f.31 (vertical dashed line). The field 
amplitude then progressively rises to its steady state value 
exp(— 2.80) « 0.06f The inset shows the similar behavior nu- 
merically obtained for aoz « 4.5 with the Brillouin parame- 
ters, namely ^/ujo = 0.071 and ujp/ujo — 1.11. 



field obtained for oqz — 2.80 with our parameters. The 
field amplitude actually cancels for jt « 1.31 and, as 
expected, the precursor consists of an unique lobe that 
clearly precedes the arrival of the main field of steady 
state amplitude e~^'^ w 0.061. The values agz « 2.80 
and jtc ~ 1.31 are obviously specific to the weak sus- 
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ceptibility and narrow resonance limit considered in the 
present paper but a numerical exploration shows that 
comparable values are obtained with theparameters con- 
sidered by Brillouin (see pp. 55-57 in [4|). A precursor 
with a single lobe is then obtained for a^z « 4.5 with 
7tc « 2.80 (see inset in Figl]). 




"1 I r 

1 2 



Retarded Time (units of 1/y) 

Figure 7: Transmitted field as a function of 7* for a^z — 
2. Inset: corresponding numerical result with the Brillouin 
parameters. 

For shorter propagation distances, the precursor, if it 
exists, is less and less distinguishable from the main field 
that becomes larger and larger. Figure [7] shows the an- 
alytical result derived from Eg. ([55]) for = 2 . There 
is only a hardly visible minimum of amplitude between 
the precursor and the main field. The numerical solution 
obtained for the same optical thickness with the Brillouin 
parameters shows a similar behavior (see inset in Figl?]). 
When aoz <C 1 (optically thin medium limit), the enve- 
lope of the transmitted pulse takes the asymptotic form 

eo{z, t)^[l- aoz (l - e"^*)] unit). (38) 

One only observes in this case an exponential fall of 
the field amplitude from 1 to its steady state value 
exp(— ao-z) ~ 1 — agz. Ea. (|38p gives a not too bad ap- 
proximation of the exact result for an optical thickness 
of 0.5 (1 for the intensity) as considered in {25| . 

It is assumed in the previous calculations that the ini- 
tial rise of the envelope of the incident field is instanta- 
neous. In real experiments, the rise time Tr is obviously 
finite and rise time effects may prevent the excitation of 
the Sommerfeld and Brillouin forerunners [19,, 20| . The 
situation is better for the precursors (22l . [251 [2g| . An ex- 
amination of Eg. ([?!]) and Eo. ([5^ in the resonant case 
indeed shows that the rise-time effects will be negligible 
if jTr -C 1 and aoZ'yTr <C 1. Only the first condition is 
usually considered in the literature. It is clearly not suf- 
ficient for the large optical thickness required to obtain 
well-developed precursors. In the millimeter-wave exper- 
iment on a molecular absorber reported in (22| . agz ~ 70 



and jTr « 3 x 10~^. We have then aozjTr « 0.2 and 
the intensity profile Jq {2^aozjt) e'^'^^unit) predicted 
by Eg. pop is well reproduced, except for the initial in- 
tensity slightly smaller than 1. On the other hand the 
optical experiments on a cloud of cold atoms reported in 
[2a | show that the amplitude (the shape) of the precursor 
may be considerably reduced (modified) when the condi- 
tion aoz^Tr <C 1 is not satisfied even when the condition 
jTr <C 1 holds. Fig. 4 in [2^ gives an example where the 
peak intensity of the precursor is reduced by about one 
order of magnitude with a second lobe larger than the 
first one. 

In the spirit of the pioneering work by Sommerfeld and 
Brillouin, we have considered up to now step-modulated 
incident pulses. The case where the step unit) is replaced 
by a single-sided exponential e~^^UH{t) is considered in 
0, nil . In the optical experiment reported in [23] , such 
pulses are obtained by passing ultra-short laser pulses 
through a Fabry-Perot resonator and the rate F is the 
damping rate of the resonator. In the nuclear forward 
scattering experiment reported in F is the decay 

rate of the 14.4 keV-state of ^^Fe used as source. In 
both cases F <§C wq and the rotating wave approximation 
holds. The envelope 60(2;, t) of the transmitted field is 
easily determined by replacing unit) by e~^'^UH{t) in the 
calculations leading to Eg. ([55]) . We get: 




X e-^*UH{t) (39) 

A very simple result is obtained in the nuclear for- 
ward scattering experiment where the source and the ab- 
sorber are made of the same material. We have then 
F = 7 and the envelope takes the form eo{z,t) 
Jo {2y^aozjt) e^^^uuit), exact whatever a^z may be. 

New experiments of resonant nuclear forward scatter- 
ing in optically thick samples were achieved in the 1990's 
by using synchrotron radiation instead of radioactive 
sources. For a review, see for example (291 . The gen- 
erated transients were named dynamical beats. Their 
theoretical study is very simple in the case of a single 
line. Indeed the duration Tp of the synchrotron pulses 
used in these experiments is long compared to I/wq but 
extremely short compared to I/7 and very short com- 
pared to 1/ [a^z^) (typically 10^ smaller). The rotating 
wave approximation is then justified and the convolution 
product of Eg. ((32)) is reduced to 

eo (.,0 « fit) aoz^A^'^'^/^%-^'uH(t). (40) 

^a.QZ^t 

Here A = JJ^ f{t)dt is the area of the incident field 
envelope (to distinguish from that of the incident field 
itself considered in SecHIl. This result is consistent with 
that given by Eq.(2.4) in [29] and with the experimental 
observations (see Fig. 2 and Fig. 3 in this reference). 
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Experiments were also achieved in optics by using 
ultra-short laser pulses. See, e.g., [2l|, [H, [5l|. In the 
experiments on a semiconductor crystal (exciton tran- 
sition) [5l| and for the largest optical thicknesses con- 
sidered in the experiments on an atomic vapor |2ll [23 | , 
the condition Tj, <C 1/ (aoz^) is not satisfied. The enve- 
lope of the transmitted pulse (polariton beat or Ott pulse) 
then differs from that given by Ea. (|40p and can only be 
determined by numerical calculations of the convolution 
product /io(z,t) (8) fit). It should however be noticed 
that, even when aozjTp = 0(1) as in the case considered 
Fig. 3c in the solution given by Eq. fiOl) fits fairly well 
the exact solution for retarded times exceeding a few Tp. 
This explains in particular why the successive minimums 
of the transient observed in [51'| at large enough retarded 
times occur at the times predicted by Eg. ipU)) . In agree- 
ment with Crisp [4l|, we emphasize that, in all cases, the 
observed transient cannot be identified to the Sommer- 
feld forerunner as imprudently stated in [52] . This erro- 
neous claim originates from confusion of the field given 
by Eq. (|14p in the limit aozj/ujQ — >■ oo (Sommerfeld fore- 
runner) with the envelope of the field given by Eg. pip 
and Eg. (15^ when aozj/ujo ^ 1. See also [ssj . 

V. CONCLUSION 

We have studied in detail what become the Sommerfeld 
and Brillouin forerunners generated by an incident step- 



modulated pulse in a single-resonance absorbing medium 
when the propagation distance decreases. Analytical cal- 
culations, combining direct Laplace-Fourier integration 
and basic saddle point method, have been made possible 
by considering the double limit where the resonance is 
narrow and the medium susceptibility is weak. We have 
shown that the structure of the transmitted field only de- 
pends on loq/j and aoz (wq resonance frequency, 7 reso- 
nance width, agz resonance optical thickness). The Som- 
merfeld and Brillouin forerunners are well apart for prop- 
agation distances z such that uqz ^ ujq/ (27^) (FigllJ, 
they overlap but remain discernible if aoz ^ ^0/7 
(Figl2]) and become practically indiscernible (Fig[3]) when 
aoz — 0(0^0/7) • Finally, they originate clean beats 
(FigHland FigE]) when 1 <^ a^z <^ wo/7. These beats 
are nothing else than the optical precursor or dynamical 
beat actually observed in various spectral domains. A 
remarkable feature is obtained for agz = 2.80 irrespec- 
tive of the value of ujq/j. The precursor then consists in 
a unique lobe clearly preceding the establishment of the 
steady-state field (Figl6]). For shorter propagation dis- 
tances, the precursor, if it exists, is less and less distin- 
guishable from the steady-state field that becomes larger 
and larger (Fig[7]). All our analytical results on optical 
precursors are obtained without making the slowly vary- 
ing approximation and are general. They are applied to 
other modulation schemes that the step modulation, in 
particular to revisit the dynamical beats, polariton beats 
and Ovr pulses generated by ultra-short incident pulses. 
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